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Abstract 



> 

■ The experimental observation of intense light emission by acoustically driven, 

periodically collapsing bubbles of air in water (sonoluminescence) has yet to 
receive an adequate explanation. One of the most intriguing ideas is that 
I the conversion of acoustic energy into photons occurs quantum mechanically, 

' through a dynamical version of the Casimir effect. We have argued elsewhere 

that in the adiabatic approximation, which should be reliable here, Casimir 
Q^. or zero-point energies cannot possibly be large enough to be relevant. (About 

D i 10 MeV of energy is released per collapse.) However, there are sufficient sub- 

' tleties involved that others have come to opposite conclusions. In particular, 

it has been suggested that bulk energy, that is, simply the naive sum of ^Huj, 
^ ■ which is proportional to the volume, could be relevant. We show that this 

^ . cannot be the case, based on general principles as well as specific calculations. 

In the process we further illuminate some of the divergence difficulties that 
plague Casimir calculations, with an example relevant to the bag model of 
hadrons. 

78.60.Mq, 42.50.Lc, 12.20.Ds, 03.70.+k 



Typeset using REVTeX 



*e-mail: milton@mail.nhn.ou.edu 
^e-mail: ng@physics.unc.edu 



1 



I. INTRODUCTION 



One of the most intriguing phenomena in physics today is sonoluminescence [^^. In 
the experiment, a small (radius ~ 10~^ cm) bubble of air or other gas is injected into water, 
and subjected to an intense acoustic field (overpressure ~ 1 atm, frequency ~ 2 x 10^ Hz). If 
the parameters are carefully chosen, the repetitively collapsing bubble emits an intense flash 
of light at minimum radius (something like a million optical photons are emitted per flash), 
yet the process is sufficiently non-catastrophic that a single bubble may continue to undergo 
collapse and emission 20,000 times a second for many minutes, if not months. Many curious 
properties have been observed, such as sensitivity to small impurities, strong temperature 
dependence, necessity of small amounts of noble gases, possible strong isotope effect, etc. 

No convincing theoretical explanation of the light-emission process has yet been put 
forward. This is certainly not for want of interesting theoretical ideas [Q. One of the 
most intriguing suggestions was put forward by Schwinger 0, based on a reanalysis of 
the Casimir effect 0. Specifically, he proposed that the Casimir effect, first considered by 
Casimir as the force between parallel conducting plates due to zero-point fiuctuations in the 
fields [0,11, be generalized to the spherical volume defined by the bubble p|-p!T|, and with 



the static boundary conditions appropriately removed. He called this, as yet, unformulated 
theory the dynamical Casimir effect. Unfortunately, although Schwinger began the general 
reformulation of the static problem in Ref. (most of which had been, unbeknownst to 
him, given earlier |T2[), he did not live to complete the program. Instead, he contented 
himself with a rather naive approximation of subtracting the zero-point energy ^ J2 of 
the medium from that of the vacuum, leading, for a spherical bubble of radius a in a medium 
with index of refraction n, to a Casimir energy proportional to the volume of the bubble: 

E,^,, = '-^fmh(,-l). (1.1) 



(27r)3 2 V n 

Of course, this is quartically divergent. If one puts in a suitable ultraviolet cutoff, one can 
indeed obtain the needed 10 MeV per fiash. On the other hand, one might have serious 
reservations about the physical meaning of such a divergent result. 

In an earlier paper, we reconsidered the Casimir effect explanation of sonoluminescence 
T3| , p!^ . We argued there that the leading term (|1.1|) was to be removed by subtracting the 



contribution the formalism would give if either medium filled all space. Doing so still left us 
with a cubically divergent Casimir energy; but we argued further that this cubic divergence 
could plausibly be removed as a contribution to the surface energy. The remaining finite 
energy, in the presumably accurate uniform asymptotic approximation, 

64^' ^^-^^ 

is at least ten orders of magnitude too small to be relevant to sonoluminescence. 

The reader might object at once that all this is in the static approximation, and the 
rapidly collapsing bubbles involved in sonoluminescence are anything but static. However, 
the time scales seems favorable for a simple adiabatic approximation to be accurate. Optical 
photons correspond to a time scale ~ 10~^^ s, while the fiash duration is ~ 10"^^ s. That 
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is, the bubble changes very httle during one period of the hght emitted. Of course, there 
may be processes here occurring on much smaller time scales, so it would be highly desir- 
able to remove this adiabatic approximation, which we hope to accomplish in a subsequent 
pubhcation. 

Eberlein has also proposed a version of the dynamical Casimir mechanism (perhaps 
more properly called the Unruh |jT6| mechanism) which she claims can explain the observed 
radiation. We have criticized her calculation on technical grounds , but mostly on the ba- 
sis of her use of ultrarelativistic velocities. See also Ref. jl^. If, in fact, reasonable numbers 
are used in her result, the energies involved are too small by 18 orders of magnitude, and 
even if her ultrarelativistic velocities are used, only lO^'^ MeV is available. So, qualitatively, 
her results are not inconsistent with ours. 

However, recently there has been a proposal that, indeed, the bulk energy result of 
Schwinger is relevant (of course, it's correct) |T^. These authors make a great issue of 
the subtraction of the uniform medium contribution, implying, it would seem, that we were 
unaware of what we were doing. Since this is a serious issue with experimental consequences, 
and since, admittedly, there are subtle issues of principle involved here, in this paper we wish 
to return to this point and provide further evidence for our result ( |1.2| ). In the following 
section we will explain more fully why this subtraction was made, indicate that it has a 
rather long history in Casimir effect calculations, and was in fact made by Schwinger in 
before he abandoned that effort. Then, in Section III, we recall the old connection between 
the Casimir effect and van der Waal forces, and show, in fact, that a finite energy of the same 
magnitude as the Casimir energy ( |1.2|) can be obtained from the latter. Finally, motivated by 
recent work on regulating Casimir energies by continuing in the number of space dimensions 
T9| , we examine, in the Appendix, whether dimensional continuation can be used to give an 



unambiguously finite value for the Casimir energy for a bubble in a dielectric, for example. 
The negative answer to the latter question shows that the quartic and cubic divergences 
found there are real. Again, appropriate physical arguments must be used to show that they 
are not relevant to the situation at hand. 



II. DEFINITION OF THE CASIMIR ENERGY 



In |T^ we derived a formula for the Casimir energy due to electromagnetic field fluc- 
tuations in a space divided into two parts by a spherical surface of radius a. The interior 
region, r < a, the inside of the bubble, has permittivity e' and permeability /i', while the 
exterior region, r > a, the outside of the bubble, has permittivity e and permeability n. 
We initially ignore dispersion. (Although it can be included [^], dispersion turns out not 
to affect our conclusions |13[-) We calculate vacuum expectation values of field products in 
terms of Green's dyadics for the corresponding classical electrodynamics problem: 



^(E(r)E(r')) = r(r,r'), 

z(B(r)B(r')) = -^V x r(r,r') x V', 



(2.1a) 
(2.1b) 



where T is the Green's dyadic for Maxwell's equations pT|. The result for the Casimir 
energy is 
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1 /-oo oo (1 

E = -— / dye'y'Y.i^l + l)x—lnSi, (2.2) 
Ana J-oo ax 

where 

= [s^{x')e[ix) - s[ix')ei{x)r - e[siix')e[{x) + s[{x')eiix)]\ (2.3) 

with 

— 1 

which is expressed in terms of modified Bessel functions 

si{x) = ^/xIl+l/2{x), (2.5a) 
ei{x) = y/xKi+i/2{x). (2.5b) 

The expression for the energy is regulated by the insertion of a Euchdean time-sphtting 
parameter, 6 = (x4 — x'^)/a, and the variables are 



It is completely manifest that (|2.2|) does not have a well-defined limit as 6 ^ — it is 



quartically divergent. Indeed, it is easy to show as does, that the quartically divergent 
term here corresponds precisely to the Schwinger result ( |1 . 1| ) when e' = /i' = 1, = 1. 
However, it is also quite clear that the calculation is not yet done when we have reached 



this point. As we stated in |]13l, "We must remove the term which would be present if either 



medium filled all space (the same was done in the case of parallel dielectrics I^TJ)." When 
we look at the latter reference, we see immediately the point. Again to quote, this time from 
pT| : "These terms [to be subtracted] correspond to the electromagnetic energy required to 



replace medium 1 by medium 2 in the displacement volume. (Since this term in the energy 
is already phenomenologically described, it must be cancelled by an appropriate contact 
term.)" What we were saying there, in the present context, is that the term in the energy 
corresponding to the boundary-condition-independent Green's function 

f['^ =zkji{kr^)h\'\kr^), (2.7) 

must be removed, because it contributes (a formally infinite amount) to the bulk energy of 
the material, which is already phenomenologically described in terms of its bulk properties. 
In fact, we are not creating material, e.g., water, we are simply displacing it when we 
insert the bubble, and force the bubble to expand and contract. The energy per unit 
element of medium is therefore not changed. (The density of the air in the bubble of course 
changes greatly, but the zero-point energy of that relatively dilute medium is certainly 
insignificant because n ^ 1. In any case, the effect of this density change is also included in 
the phenomenological description.) 

Indeed, the spectacular agreement between the the Lifshitz theory of parallel dielectrics 
1 2^, rederived in [^, and the beautiful experiment of Sabisky and Anderson p3[ seems 



strong vindication of this subtraction procedure. 
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Further evidence that we are on the right track is provided by Schwinger himself. In 
the first National Academy article cited in [§, where he rederives the result for parallel 
dielectrics, he explicitly removes volume and surface energies: 

one finds contributions to E that, for example, are proportional ... to the volume 
enclosed between the slabs. The implied constant energy density — independent 
of the separation of the slabs — violates the normalization of the vacuum energy 
density to zero. Accordingly, the additive constant has a piece that maintains the 
vacuum energy normalization. There is also a contribution to E that is propor- 
tional to [the area], energy associated with individual slabs. The normalization 
to zero of the energy for an isolated slab is maintained by another part of the 
additive constant. 

Admittedly, the situation is more clear-cut in the parallel-plate geometry. However, in the 
following paper (the last reference in ^) where Schwinger begins to set up the problem 
for the spherical geometry (but leaves the details to Harold ||2^), a close reading shows a 
similar subtraction is implicit. Unfortunately, when Schwinger went on to apply Casimir 
energy to sonoluminescence in 0, he does not make use of the general analysis in 0. Our 
interpretation is that at that point Schwinger lost the energy or courage to complete the 
full calculation, and needing an immediate result to confront the phenomenology, simply 
jumped to the unsubtracted, unregulated result (|1.1| ) — see the second reference in 



But enough of argumentation. Let us turn to detailed calculations that support our 
contention. 



III. DERIVATION OF CASIMIR EFFECT FROM VAN DER WAALS FORCES 



It is familiar that the van der Waals forces between polarizable molecules — the Casimir- 
Polder forces [^] — can be derived from the Casimir forces between dielectric bodies. We 
interpret this as meaning that the Casimir effect is merely a local field form of the action- 
at-a-distance summation of the forces between the molecules that make up the material 
bodies. 



Let us begin with a variation of the argument given in ||2T|. Consider a dielectric slab 
bounded by planes z = and z = a, having dielectric constant e; outside this region there 
is vacuum, e = 1. According to the Lifshitz formula 121,22|, the force/area between the 
surfaces is 



/ = - 



2tt Jo 2tt 



\K3- Ki 



- 1 



-1 



(3.1) 



where, in the ith medium (we denote the region of the slab by 3, that the outside regions 
by 1), 



K 



(3.2) 



Now suppose the medium is tenuous, so that the dielectric constant differs only slightly from 
unity. 
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e-Kl. (3.3) 
Then, with a simple change of variable, 

= (p, (3.4) 
we can recast the Lifshitz formula (13.11) into the form 



^""3^ rd(('rp(0-irW-ir + l]e-'''^. (3.5) 

If the separation of the surfaces is large compared to the characteristic wavelength charac- 
terizing e, a(c ^ 1, we can disregard the frequency dependence of the dielectric constant, 
and we find 

^ _23(.- iP 

6407r2a4 ^ ^ 

For short distances, a(c ^ 1, the approximation is 

1 1 fOO 



327r2 a3 Jo 

These formulas are identical with the well-known forces found for the complementary geom- 



etry in 21 



Now we wish to derive these results from the sum of van der Waals forces, derivable from 
a potential of the form 

y=-^- (3-8) 

We do this by computing the energy (A^ =density of molecules) 

E = -iB/V'^°d^/;d.'/(drJ(*'J^^— (3.9) 

If we disregard the infinite self-interaction terms (see below), we get 

dE_ 1 
^ da A (2-7)(3-7)a7~3' ^ ' ' 

So then, upon comparison with ( |3.6| ), we set 7 = 7 and in terms of the polarizability, 
we find 

B = f a^ (3.12) 

or, equivalently, we recover the retarded dispersion potential. 
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V--'^%, (3.13) 



whereas for short distances we recover the London potential, 



Q 1 /-oo 

y=—J dCaiO'. (3.14) 



IT r" Jo 



Our intention is to carry out the same simple calculation for a dielectric sphere. The 
first couple of steps are unambiguous (6* is the angle between r and r'): 
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E = -^BN^ l{dr){dr') 

a r 1 1 1 

(3.15) 



2-7 Jo 



dr / dr' rr 



^2 _^ ^12 _ 2rr' cos O)"!/"^ 

1 1 



Now, however, there are divergences of two types, "volume" (r' r) and "surface" [r —>■ a). 
The former is of a universal character. If we regulate it by a naive point separation, r' r+6, 
(5 — > 0, we find the most divergent part to be 

ttBN^ 1 Ana^ 

which is identical to the corresponding (omitted) divergent term in the parallel dielectric 
calculation, where V = aA. This is obviously the self-energy divergence that would be 
present if the medium filled all space, and makes no reference to the interface, and is therefore 
quite unobservable. This is the analogue (although the e dependence is different) of the 
volume divergence in the Casimir effect, ( p,.l| ). 

If, once again, the divergent terms are simply omitted, as may be weakly justified by 
continuing in the exponent 7 from 7 < 3, we obtain a positive energy, 

23 

E^^y, = -——{e-l)'. (3.17) 
153o7ra 

This may be more rigorously justified by continuing in dimension, a procedure which has 
proved useful and illuminating in Casimir calculations |]Ty|. Thus we replace the previous 
expression for the energy by 

E = --BN^ I d'^rd'^r'- — (3.18) 
2 J |r-r'|^ ^ ' 

where, in terms of the last angle in D-dimensional polar coordinates. 



/ 



d^r = /„ / t^^sin^-2 0. (3.19) 



r 



If we take, say, r' to lie along the z axis, so that 9 is again the angle between r and r', we 
find 
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E = — BN 



dr r 



ID-I 



drr 



D-l 



2 r (f ) r (-^) Jo 

X c/cos^(l-cos2^)(^-=^)/'(r2 + r'2-2rr'cos^)-^/2. 



(3.20) 



The angular integration can be given in terms of an associated Legrendre function P^[z), 



dt (1 - t2)(D-3)/2(^2 ^ ^,2 _ 2rr't)-in 



^2 I ^'2 
I ^ 2 ^_ ^ ,^/2 1 



(3.21) 



Now let us substitute this into the expression for the energy, and change variables from r, 
r' to 



2 /2 

z = r + r , y 



y^2 I ^'2 

Ly^ 2 ^/2 



(3.22) 



The X integral is then trivially done, leaving us with 



E 



'2 Ji 



dy 



2a 



2 \ -0-7/2 



{y'-ir-''"Pl;%2{y\ (3.23) 



valid for D > 7/2. Integrals of this type are given in ||26|| : 

/CO 
dy {y - ly'^/^y + if+'^/'-'P.^iy) = 2^— - 



T{-2b) 



(3.24) 



6-a)r(l-6)' 

valid for Re a < 1, Re 6 < 0. Then we have, using the duplication formula for the F function. 



E = -BN' 



r(D/2)r(D-7/2 + i)(z^-7)' 



(3.25) 



The resulting formula is regular when D and 7 are both odd integers, so we can analytically 
continue from D > 7 to 2? = 3 for 7 = 7. Doing so gives us, using Eq. ( ^.12| ), 



E = BN' 



.TT^l 23(e-ll 



24 a 24 647ra 



(3.26) 



exactly the same as ( |3.17| ). Note that the magnitude of this result is nearly the same as that 
found in [|l^], and stated in (pT^) , differing only by the factor 



23 4 

= 1.22, 

24 TT 



(3.27) 



which is a plausible difference in that the previous calculation was only in the leading 
asymptotic approximation, but the sign is opposite! We offer as evidence for the validity 
of this methodology the fact that the formula ( |3.25| ) gives the correct Coulomb energy of a 
uniform ball of charge, for which 7 = 1. 
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Evidently, we have reached the frontier of our understanding of the Casimir effect and its 
connection with van der Waals forces. The subtraction procedure may well be ambiguous, 
although the volume and surface divergences are unambiguous. That these divergences 
are real is further reinforced by the considerations of the Appendix, which shows that the 
technique of dimensional continuation fails for this case. But these divergences are not 
relevant to the light emission process, although they would be to a first-principles calculation 
of the energy density and surface tension of the medium [^|. However, our qualitative 
conclusion, that quantum vacuum energies are completely irrelevant to sonoluminescence, 
is dependent only on the order of magnitude of the finite remainder, given by either ( p..2| ) 
or (CT - 



IV. CONCLUSIONS 

Our conclusions here are threefold: 

• The divergences that occur when interior and exterior modes are mismatched, whether 
by exclusion of one set, or by changing the the speed of light in the two media, are 
real, and cannot be circumvented by a mathematical trick. 

• Volume divergences are not physically meaningful, since they reflect self-energy effects, 
and serve to define the intrinsic properties of the material. They are naturally can- 
celled out by the introduction of a suitable contact term. What is left is a surface 
divergence, which presumably is physically meaningful, yet should be absorbed into a 
renormalization of physical parameters, such as the surface tension. 

• The magnitude of the finite remainder, of order 1/a, apparently may be extracted un- 
ambiguously. Whatever its sign, it is far too small to be relevant to sonoluminescence. 



ACKNOWLEDGEMENTS 



We thank M. Visser for discussion about his work and the significance of the bulk energy 
contribution. We thank the U.S. Department of Energy for financial support of this work. 



APPENDIX: DIMENSIONAL CONTINUATION OF THE CASIMIR EFFECT 

The fact that the above dimensional regularization of the van der Waals energy gave 
a finite result suggests that we re-examine the Casimir calculation to see if possibly an 
unambiguous finite result could thereby be obtained. We will not be surprised to find a 
negative answer to this question, since the perfect cancellation between interior and exterior 



modes cannot hold true with different speeds of the light in the two media p7 |. 

We will content ourselves by examining the extreme case of e — >■ oo in the exterior region, 
that is, a bag with perfectly conducting boundary conditions on the surface. Since it is 
necessary to continue the individual modes, we will examine the TE mode as representative. 
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[As we will see, the subleading divergences cancel between the TE and TM modes.] In three 

x^^, (Al) 



dimensions the interior modes alone give |28 



n=l 



SnXX) 



where the generalized modified Ricatti-Bessel functions are 



Sn[X) 



X 



x) 



X 



(A2) 



where v = n + D/2 — 1 (=72 + 1/2 here). The generalization of this result to D space 



dimensions is 19 



E: 



TE 



1 



1 



where the weight function is 



win, D] 



— Y,win,D) 

^) n=l 



2uT{n + D -2) 



dxx 



Sn{x)' 



(A3) 



(A4) 



Again, as elsewhere, in the "vacuum energy" term was subtracted. (As noted in the 
text, the justification was only partly that it removed the most divergent terms.) This was 
obtained from the free Green's function, which in D space dimensions is 



^°(r,r',i 



2uT ( 



D 



fCf/2-^)(cos^)J.(fcr<)i7«(fcr>), 



(A5) 



in terms of the ultraspherical or Gegenbauer polynomial. The stress on the sphere is obtained 
by applying the appropriate differential operator corresponding to the stress tensor. 



Trr = -(V,V,. + CU^ - • Vx')e^° 



n{n + D- 2y 



(A6) 



Subtracting this from the previous result gives 



TE 



r> w(n,D) dxx{ [ 



+ X 



3-D 



s„(x)e„(x) 



1 + 



n{n + D-2y 



(A7) 



The question now is whether the continuation procedure described in can be success- 
fully applied here. There, we first made the integrals convergent by adding a suitable term 
to the summand which sums to zero for sufficiently small dimension. Here this suggests that 
in the above integral we replace 



10 



Sn{x) 



ax ax 



(A8) 



for then the large x behavior of this term 1 + (4z/^ — 1)/ (8x^) + . . .. The vacuum subtraction 
term cancels the leading term here, leaving for the leading term in the braces in 



iD-3){D-2) 



+ 0(x- 



(A9) 



So, not surprisingly, the integral is still, in general, logarithmically divergent, although for 
D = 3 or 2 it does converge. 

Therefore, it appears that we cannot meaningfully continue off the integers. So we are 
forced to retreat back to D = 3. There we have 



D=3,TE 



where Q„ is the convergent integral. 



27ia 



^(2n + l)Q„, 



(AlO) 



n=l 



roo ( ^ ^ 

Qn = dxx < — \n \/2TixIy{x) + c.t. > . 
Jq \dx J 



(All) 



If we use the uniform asymptotic expansions for the Bessel functions, we easily find 

^3(2^2 + 3)- 



Qr 



dz z 



ivz 



16 



oo 



(A12) 



where t = (1 + z^)~^/^. The leading term here is precisely the negative of that found in the 
TM mode; that is, for electrodynamics, or linearized QCD, the interior modes contribute 
a quadratically divergent sum, rather than the cubically divergent one due to each mode. 
Practical methods of dealing with this divergent Casimir energy, which is relevant in hadronic 



physics, were suggested in |29 
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